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NATIONAL ADVISORY COmilTTEE FOR AERONAUTICS. 
TECHinCAL MEMORAl^DUM NO. 493. 

aerodyiia::ic characteristics of thin empirical profiles AIID ■ 

THEIR application TO THE TAIL SURFACES AND AILERONS OF AIRPLAl^SS.* 

By A. Tous saint and E. Carafoli. 

The increasing use of airplane wings equipped with ailer- 
ons and the importance of knowing the aerodynamic chaxacteris- 
•tics of tall units (empennages) provided with movable parts 
(rudders anc elevators) impart some interest to the so-called 
"empirical profiles." This term is applied to wing profiles 
which 3-re not drawn according to any theoretical method. Gener- 
ally these profiles are .of no particular interest, since they 
can al^veys be replaced by theoretical profiles which satisfy 
the sa;vte aerodynamic and structural conditions. For the above- 
mentioned applications, however, the profiles, modified by the 
deflection of the aileron, elevator or rudder, necessarily be- 
come empirical profiles, of ¥/hich it is important to know the 
aerodynaumic characteristics and, above all, to know how these 
characteristics are- affected by the magnitude of the deflec- 
tions. 

For this purpose we thought best to employ the method pro- 
posed by ilunk for the approximate theoretical study of thin, 
slightly curved profiles, assimilable, f rom an aerodynamic view- 
*From L'Aerophile, June, 1928, pp. 179-183. 
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point, to their mean cai^iber line. 

Principle of the Method 

Given a profile isrhose mean camber line ASB is moder3.tely 
curved and v/hose chord AB = I. . 

It is v;ell known that, in the conformal transformation 

Z = S + ^ (1). 
I 

a circle C, of radius ^ - 9-^^ centered at 0, gives a 
straight line AB of length I = 4a. 

On a,dopting the above transformation function (l), the 
line ASB, very close to the chord AE, will have for its 
antecedent in the ^ plane a curve S' which will lie very, 
close to the circle C, eaid which will cut this circle at the 
points A' rjid 3' , antecedents of A and B. 

The curve S' rnay be defined by the expression 

r = a (1 + r) e^^ (2) 

in which r is the radial segment, as a function of G, meas- 
uring the sirall distance betweai the circle C and the curve 
S' along the radius vector of " a,rguiiient " S. 

The mean line ASB will be obtained by the transformation 

z = e' + (!') 

which, on neglecting the terms in r^, gives 
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Z = 2 a COS0 + 2 1 ar sin 6 . 
For the abscissa X = 2 a cos 6 

} (3) 

the ordinate vrill be Y = 2 a r sin0 

♦ 

To transform the circle e into the curve S' so that a 
point t,' = a (l + r) e^^ of the curve S' will correspond to 
a point S = a e'^'^ of the circle 0, let 

= e + p (4) 

Under these conditions, the trailing edge B of the pro- 
file ASB will correspond to the point B' , (9 = tt) of the 
curve S' and to the point B'q (cp = tt + of the circle C. 

According to Joukov/ski's hypothesis, the point B'o must be a 
point of zero velocity. Hence, the circulation will have the 
value 

T = 4 TT a V sin (a + P^) (5) 
The transf orination function under the general form 

r = S + a (Ao + i Bo) + a? ^ + • . . 

+ a 5 (6) 

^n 

will enable us to pass from the circle 0 to the curve S'. 
On replacing t,' and t, by their respective values, 

^' = a (1 + r) e^® and ^ = a e"^"^ 
and, since P and r are very small, 
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r ^ S An' cos (n + 1) e + S Bn sin (n + l) 9 
P = E An sin (n + l) 6-2 Bn cos (n + l)e 

Since the profile ASB has no thickness, the ordinate 

y = 2 ar sin 6 

takes the same value for ± 9 . Hence, r changes sign with 9 
and 

Ao = Ai = Ag = = An.- 0. 

This gives 

r = S Bn sin (n + 1) 8 ^ 

} (7) 

3 = - Z Bn cos (n 4- 1) 6 

The coefficients Bn are therefore defined by the formula 

Bn = / r sin (n + l) 9 d 9 (8) 

The aerodynajnio characteristics of the profile ASB are then 
given by the following formulas: 

1, Lift P = p V r = 4 rr a p sin (a +p) 

whence = 2 tt (a + pg) . (9) 

2 . Moment = Mq - P 2 a co s a , 

Let Mj^ = 4 rr p T= \%^\ . (a + Y ) - 8 tt a^ p (a + ^o) 
whence G^.^^>^ ^ | (a + 7 ) _ tt (a + Po) (10) 

Xi =0^ l^"^*^ being the term in i of the transformation func- 
tion (6). Thus, we have 
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Xi = e^^'^ = (1 + i Bi) (11) 

The computation of and. '^m(_^) leads therefore to the 

calculation of Pq, \'LA and of "Y . 



Calculation of 

This angle corresponds to the point Bq^ , that is, to 
e = tr . Under these conditions the formulas (7) give 

3o = - 2 Bn cos (n + 1) TT = Bq - B^ + Bg - Bg + ... 

Replacing Bq , B^^, Bg ... etc., by their values derived from, 
formula (8), we obtain 

= / r (sin e - sin 26 + sin 3 9 ...) dO 



. ie 

^ e 1 + e^^ 

which reduces to 

P° = TT-^e 21 + cos e 2-TT-i: ^0 y 1 + cos e ^^^^ 

In practice the profile abscissas and ordinates are re- 
ferred to. the axes x' A y' and the relative coordinates 



_ y 

4a 

are preferably considered. 



~ 4a " 4a 
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Thus equation (12) gives for the value 



'^ =i/^— — — = Ti f, {%) d i (13) 
° (1 -i)Ji (1 e) ° 



by letting 



(e)= i (I3t) 

TT (1 - i) J% (1 -1) 



When the mean line ASB of the profile is known, equation 
(13) cam be integrated "by employing computed values of f.^ (0* 

Calculation of (Xi) and of 7 

Since is small, formula (11) gives 

= |zj = a3 (1 + B^^) S a^ 

Similarly tan 2 7 s 2 7 = Bi 

The value of Bj_ is found by equation (8), which gives 

-1 sir o TT 

Bi=^/ rsin2ede=-£./ ycosGde 

Letting = ^ and employing the relative coordinates 
^ and r|, we obtain 



^J-o = 2/ -n (1 - 2 e) = 2 /' Tif^ (e) d I (14) 

y I (1 - e) 

with fc (n '= - t " " ^ (14') 

y e(i -e) 
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■ As before, the values of f g ( ^ can "be calculated for 
different' abscissas I,- aiid' th^ ' integral; ion of equation (14) 
will be accomplished according to the values of r\ obtained on 



the mean camber line ASB. ' Expression (10) for ^m(^^'^ then 



becomes . 



-'m 



(A) , 2 \ 2 7 



- TT (a + 3o) = - 0.25 O2 + 0 



- R R 

- A ^1 " S PO 



(10«) 
(10") 



In short, the aerodynamic characteristics of empirical 
profiles, assimilable to their mean camber line, due to their 
relative thickness and curvature, can be calculated from the 

formul8.s 

C2 S 2 TT (pc + Pq) 2TT [a + T) f, (e) d e] (9') 



% - 0.25 O2 + Om^ = - 0.25 



+ 2 / T] fg (1) d e -2 / Ti f^ (I) d e 

o 0 



(10'") 



for which the function f^ (i) and fg {i) are given in the 
table beloTif e.s functions of i . 



e 




0.025 


0.05 


0.10 


0.20 


0.30 


0.40 


0.50 




it) 


2.09 


1.54 


1.18 


1.00 


0.99 


1.08 


1.27 




it) 


6.10 


4.13 


2.67 


1.50 


0.87 


0.41 


0 , 


i 




0.60 


0.70 


0.80 


0.90 


0.95 








(1) 


1.62 


2.31 


5.98 


10.6 


29.2 
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ii) 


-0.41 


-0.87 . 


-1.5 


- 2.67 


- 4.13 







IT.A.O.A. Technical Memora.nd.um No. 493 



8 



Profiles Having Any HeaJi Camber Line 

If the nea.n camber line of the profile cannot be represent- 
ed by a simple relation between "n and i, the integrations 

fi (o d t s-nd /' Ti (e) cL e 

0 (5 

can be ma,de either graphically or mechanically, by making the 
products T] and Ti fg according to the values of r| derived 

from the niean camber line. In general, the value of T) f^^ in 
this calculation increases indefinitely toward the trailing 
edge (i.e., for ^ = l). We can then evaluate the integral 
n fj d C , from 1=0 to | = 0.95, and estimate the por- 

0 

tion of the integral from | = 0,95 to | = 1.0. Assuming 
that this portion of the mean camber line is rectilinear, v/e 
find that the additional contribution amounts to 2.9 r\', 
being the relative ordinate for | = 0,95 (according to Glauert). 

Profiles with Particular Mean Camber Lines 

Fnen the form of the mean camber line can be expressed by 
a simple ratio between r\ and I , the values of ^''o 
can be calculated by direct integration. For example, the 
equation 

T] = b e (1 - e) (c - e) 

represents a profile whose mean camber line is of double curva- 
ture when c varies between ^ and 1.. Direct integration gives 
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Po = I (4 c - 3) 

- 32" ' • 

from which « _ rrb 



= If (7 - 8c) is derived. 

Application to Taal Units with Movable Parts 



For c = ^ we get Om_ = 0. 



For a stabilizer having a symmetrical biconvex profile, 

the deflection of the elevator HB by a certain moderate ajigle 
P, forms a profile whose mean camber line is composed of two 
straight lines, AH aJid HB. 

Let a be the ra.tlo of the movable part HB to the total 
surface area AH + HB, then: 

a - HIB. e, , ~ Mi - 1 „ a and h = ^ = (l - a) B 
A B' (H) =^ 3 - J- - " ana n ^ g \i. p 

For I varying from 0 to (l - a) we get 

and for | varying from (l - a) to 1 we get 

ri=|(l-e). . ,. 

Under these conditions the values of Pq and IJ-q are 
calculable by direct integration of formulas (13) and (14). 
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Thus we get 



TTo 1 - a (1 _ 1) 71 (1 -i) 



(1 - i)J% (1 ~ i) 



1 - ^ /I (1 - e) 



and 



and finally 



^2/ , h (1 - 3 e) di 
1-0 e(i~i) 



2h 1 _^ J]_ arc cos//a 

^ |y5-(i - a) 2a - a (1 - a)_ 



(15) 



, cr , f 2a - 1 , TT arc COS\/a l /'^c;^ 

and = h 1 + _ jy (16; 

j^ya (1 - a) 2a a (l - a) 

Aerodynamic Characteristics 

Here we have: 

O2 = 2TT (a + fc). 

In practice, the lift coefficient is expressed as a func- 
tion of the- deflection P of the elevator and of the angle of 
attack of the relative wind on the stabilizer. AH. Therefore, 

we may write 

Oz = 2rr (-.a^ + mo P ) (17) 
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Since hSa (l - a), and p and a =, a^. +...\., =..Ai.. .+ ^^P, 
we finally obtain 

0,7 = 2tt [tti + fl + ^ j'o (1 - a) - ^ arc cos^/'a^^] 

^ \ TT . TT / . 

Hence 

mn = 1 - - C arc co s ^a" - {1 - a ) ] ■ ( 18 ) 

Figure 3 represents the variation of with respect to 

a according to formula (18). The experimental results of dif- 
ferent aerodynamic laboratories, for deflections P varying be- 
tween ±5° and ±15°, are here shown. The locations of these 
representative points with reference to the theoretical curve 
show that the value of the coefficient m, applicable in the 
formula 

O2 = 2tt ( a + m P ) , 

diminishes somewhat with the magnitude of the deflection 3 . 
As regards this analysis, we need only to call attention to the 
fact that the experimental points very close to the theoretical 
curve apply to elevators extending over the whole tail span 
without any central cut-out. When there is such a cut-out, the 
experimental values of m are somewhat smaller than the theo- 
retical values, and the empirical -expression, 

m==j'V (19) 
very closely represents the mean values applicable to moderate 
deflections (P < 10°). 
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Similarly the moment coefficient Cm 



is finally ob- 



tained as 




= - 0.25 0 



^ - 2 /a (1 - of p 



(20) 



The experimental data available for verifying this theo- 
retical formula are, however, very few* 



The preceding results may be extended to the case of sup- 
porting wings equipped with ailerons for the purpose of in- 
creasing the maximum lift coefficient. 

. Let ASB represent the mean camber line of a wing profile 
the rear portion of which HB forms an aileron rotating about 



the hinge H (Fig. 4). The aerodynamic characteristics of the 
undisturbed mean camber line are given by formulas (9') smd 
(lO"), as previously demonstrated. 

Let 3 represent the angle of deflection, such that the 
mean camber line becomes AHB. As previously shown, the approx 
imate theory is applicable to modera.te values of deflection P , 
so that the chord AB' of the modified profile is very nearly 
equal to the chord AB of the original profile. 

The aerodynamic characteristics of the profile ASB' will 
be calculable by the application of formulas (l3) and (14) to 
the equation of the line ASB' with respect to the axes ox' 
and oy» , 7 being the angle which ox makes with ox'. Let 



Application to Wings Provided with Ailerons 



N.A.O.A. Technical Meiaorand-um No. 493 13 

a = il tlien 7 ^ Qf P . For the portion AH, characterized by 

AB 

the, relative .abscissa ■ i = (l - a) we then get: 
^ C = C cos 7 - T] sin 7 s J 
' Ti' = T] cos 7 -!- ^ sin 7 S T) + e a 3 

In an analogous manner for the portion HB' ■ we get: 
V = i 

Ti' =i T] + (1 - a) p (1 - O 

The calculation of ' ^o^ "^^^^ "^^^^ ^'^ 

formulas (13) and (14) applied between 0 and (l - a) for 

the portion AH and between (l -o) and 1 for the portion 
HB' . We thus find 




+ (1 - a)3 r (1 - n fi (n <il (13«) 
and = 2 ^ f3 (U d. e + 2 [a P/^^ e fg (1) dl 

0 0 

+ (1 - a) P (1 - e) fg (E) d e] (14') 
1-a 

The first integral in Pq' represents the angle of zero 
lift of "the profile ASB. The other two integrals repre- 

sent the angle of zero lift of the deflection 3 of the eleva- 
tor. This sjigle being A 3o, ""'^e h3.-ve 

^zrr,^ = 2 TT (a + + A Po) 
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Po is calculable from the equation r\ = f (i) of the initial 
mean camber . line''' ASB. ' ' ' 

A Pq can be calculated as already described for the mova- 
ble eleva.tor, i.e., by the formula 

A Po = (^o - ^ 
with fflo = 1 - |. [arc cos y'o"- J^a (l - a) ] 

If the chord AB of the original profile is taken as the 
reference line, the angle of attack is then aj^=a_7^cc_op 
and the result is 

Oz = 2 n (a^ + Po + Ki 3). 

The first integral in the expression for la^' represents 
the value of P-q for the original profile. The other two inte- 
grals represent the correction A due to the deflection of 
the aileron, A |J-^ will be calcxilated as in the case of the 
elevator, and we have 

Omp(P) = ^mo(p=,o) - (5 ^o ~ ^ ^^o) 

or 

°^o(P) = ^«^oO=o) - 2 /^(l - af p 

Experiments with wings equipped with ailerons verify quite 
well the theoretical values of A Pq, in particular for posi- 
tive moderate deflections P . In general, the values of 
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under the same conditions, are somewhat larger than those found 
experimentally. The deviation is similar to the one found im 
the case of wing profiles having an appreciable Ckiq. 



Translation ■by- 
Nat ional Advisory Committee 
for Aeronautics. 
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